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The dynamics governing the evolution ol a many body system is constrained by a nonabelian 
local symmetry. We obtain a general form of the global macroscopic conditions assuring that 
at the microscopic level the evolution respects the overall symmetry constraint. These con- 
^ , ditions does not depend on detailed dynamics of the system. It is shown that a nonabelian 

internal symmetry leads to the new subsidiary condition imposed on the system. A structure 
of constraints is analysed. 
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The statistical mechanics of finite systems is a powerful tool to describe high-energy strong interaction phe- 
^ ' nomena. There is a well established tradition of this approach, starting from classical ideas of Koppe [Q and 

I Fermi |^ . To get a better insight into properties of systems under considerations even the very refined relativis- 

. tic statistical thermodynamics Q must be supplemented with a kinetic theory [^-Q . An available phase space is 

' always constrained by internal symmetries of the system. A consistent formulation of thermodynamics for with 

internal symmetries taken into account ||^,^ allowed for suitable modifications of the equilibrium distributions 
|-|ll|. In the case of the kinetic theory the situation is more involved. The influence of the internal symmetry 
on the generalized Vlasov - Boltzmann kinetic equations was considered recently by the authors ||l^ . 

In this paper we consider a more general scheme which allows to obtain a set of conditions due to the 
internal symmetry preserved by the underlying local interaction. We obtain a generic form for these conditions, 
independent on details of the kinetic equations. 
' Let us consider a multiparticle quantum system with the local interaction invariant with respect to the internal 

^ , symmetry group G. We shall call this underlying symmetry of microscopic interactions a microsymmctry of 

C ' the system. The system transforms under a given representation of the symmetry group. We shall call this 

property a macrosymmetry of the system. Let us try to answer a question: is macrosymmetry preserved during 
, , ■ a time evolution of the system? 

r> I The system consists of particles belonging to multiplets aj of the symmetry group. Particles from the given 

■ multiplet aj are characterized by quantum numbers vj — related to the symmetry group, and quantum numbers 

characterizing different multiplets of the same representation aj. The number of particles species is: 

3 Cj 

where d{a) is the dimension of the representation a. 

The number of particles of the specie {a,iya',C} is denoted here by A^^"-\(;-)- These occupation numbers are 
time dependent until the system reaches the chemical equilibrium. However the representation of the symmetry 
group for the system remains constant in the course of a time evolution. 



N!):i\---,Nt:^). All 



Let us consider a multiparticle state vector in occupation number representation: 

other variables, related to phase-space properties of the system are suppressed here. This vector describes 
symmetry properties of our systems and transforms as a direct product representation of the symmetry group 
G. This representation is of the form: 

af°^'®af°^'®..-0af°"'. (2) 
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A multiplicity A^("j ) of the representation aj in this product is equal to a number of particles which transform 
under this representation: 

^^-^-EfE<l.)l-E<;'- (3) 

j \ / 3 

The representation given by Eq. (H) can be decomposed into direct sum of irreducible representations A^,. Cor- 
responding states are denoted as |Afe, Aa^.; A/") where Aa^. is an index numbering members of the representation 
A and A/" is a total number of particles 

A/" = E ^^TJ ■ (4) 

k 

Each physical state can be decomposed into irreducible representation base states |Afc, Aa^ ; A/"; ). Variables 
^A are degeneracy parameters required for the full description of a state in the " symmetry space" . Let us 
consider a projection operator on the subspace spanned by all states transforming under representation A. 



7?^ 



This operator has the generic form (see e.g. 13|): 

P^ = d{A) [ df,{g)x^'^Hg)U{9). (6) 



G 



Here x'^^^ is the character of the representation A, d^{g) is the invariant Haar measure on the group, and U{g) 
is an operator transforming a state under consideration. We will use the matrix representation: 

C/(g)|<"^^),...,<";)) (7) 



(qi) n("i) n^""^ n^""^ 



f^'v^.v^ is a matrix elements of the group element g corresponding to the representation a. Notation convention 
in Eq. arises since there are nI"^,'^ states transforming under representation aj and having quantum numbers 
of the z^Q,^-th member of a given multiplet. 

Let P^'\ti) (c«„) denotes the probability that nI)21\ ■ ■ ■ , A,^"^' particles transforming under the symmetry 

group representations ai, . . . ,q;„ combine into Af particle state transforming under representation A of the 
symmetry group. This is equal to a norm of the vector and it is given by 

(7V(-), • • • , Ar(-) I P"^ . . . , A(-)) = J2 |C{<.>,..., <")}(^^)l' ■ 

Left hand side of this equation can be calculated directly from Eqs.@ and (|^). One gets finally 

^ " G 

where A^''^^ is a permutation normalization factor. For particles of the kind {a, C} we included in Eq. the 
permutation factor: 

The permutation factor A^^^ is a product of all "partial" factors 
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Because of macrosymmetry all weights in Eq. (p|) should be constant in time. So we have 

= 0. 



'} 



(11) 



(12) 



This condition assures that in a dynamical evolution the macrosymmetry of the system is preserved. By mean 
of Eq. the macrosymmetry conservation condition ( p^ is converted into a general time evolution condition 
of the form 



= 



dt 



(13) 



EE 



dN^'' 



dt 



All integrals which appear in Eq. (|^) and Eq. (^3|) can be expressed explicitly in an analytic form for any compact 
symmetry group. 

To write an expression for the time derivative of the normalization factor A^''^^ we perform analytic continua- 
tion from integer to continuous values of variables nI)"'^\ Thus we replace all factorials by the gamma functions 
of corresponding arguments. We encounter here also the digamma function -0 

d logr(x) 



■0(2;) 



dx 



This allows to write: 

dA^ 
dt 



1 Cj 



dN. 



dt 



(0 



dN^'''\,, 
dt 



(14) 



(15) 



These subsidiary conditions (13) and (15) are the necessary conditions for preserving the internal symmetry on 
the macroscopic level. Rates of change dN'^^.f^^^/dt are related to "macrocurrents" , which are counterparts of 
"microcurrents" related directly to a symmetry on a microscopic level via the Noether theorem. Eqs. ([l^ ) and 
( p^ ) can be considered as a set of conditions on macrocurrents to provide consistency with the overall symmetry 
of the system. 

For multiplicities encountered in high energy heavy ion collisions processes all digamma functions in Eq. ( p^ ) 
can be replaced by corresponding logarithmic functions according to the asymptotic formula Jlif : 



V'(iV + 1) ~ \ogN 



This gives a simplified form of Eq. ( |15| 



dt 



-^^^^EE 



E 



dN. 



(16) 



(17) 



One obtains then from Eq.(|l^) 

n 

EE 



^^(.,)diog^;^t,. 



dt 



dm: 



(18) 



i 



where 
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^ " G 

is analytically extended for continuous values of variables N^"' ■* . 

Eqs ( p^ and (18) are meaningful only for nonzero values of coefficients (p^. Let our symmetry group be a 
Lie group of rank k. It means physically that there are exactly k simultaneously measurable charges related to 
the group G. It is easy to see that weights (^) are nonzero only for such values of -/vi"^"* which are consistent 
with conservation of those charges. For the isospin SU{2) group that is the third component of the isospin, 
for the flavour SU (3) that would be the third component of the isospin and the hypercharge. In general case 

(a ■) 

one has k linear relations between variables Ni,J, what reduces correspondingly the number of independent 
variables. 

New results are obtained only for nonabelian symmetries. The local abelian charge conservation on the 
microscopic level is equivalent to the charge conservation of the multiparticle system. In the case of the abelian 
symmetry group all representations are one-dimensional ones, labelled by the charge value a. All related 
permutation factors (^0|) are equal to 1 and Eq.(||) gives 





We see that for abelian internal symmetries Eq.(p^) gives an obvious result 

For nonabelian symmetries one gets similar relations for charges belonging to the center of the group G. 

The "nonabelian" condition ( p^ ) is much more involved. The time derivatives dN^'^,^^-^/ dt can be obtained 

from analysis of the microscopic basic interaction. The most general approach ||l^ goes through Heisenberg 
equation 

'-^=AhM (22) 

averaged over states NvZl \ ■ ■ ■ ' ^'^a'^^'^ ■ ^ denotes here a particle number operator. 

Another way is to use semiclassical kinetic equations for distribution functions /^^"^"'^'■'(r, r, t). The variables 
(r, r) denote a set of phase - space variables. The number of particles of the specie {a, Va] Q 



= j di/dr47°'(r,f,t). (23) 



The case of the generalized Vlasov - Boltzmann kinetic equations was considered in 



;(C) 



In general the derivatives d N^°'\^s/dt are expressed as functionals over variables 



dt --^--(0 



(24) 



Substituting this into Eq. ( p^ or Eq. ( p^ ) one gets a new relation between variables nI'^'\ This leads 
in principle to the further reduction of the number of independent variables. Because of a highly nonlinear 
character of the resulting constraint it would be too optimistic to expect an analytic solution of this constraint. 



However, it is possible to get at least numerical estimations 15 1 based on kinetic equations ( |24| ) and taking into 



account k solvable constraints related to the center of the group G and the implicit constraint resulting from 



Eq. (|13|) or Eq. (|18|). 

There is a simple geometrical interpretation of our problem. Solutions of Eqs (|2j) span M^-dimensional 
manifold parameterized by the time variable, with W given by Eq. (|^). Constraints reduce this manifold to the 
{W — k — l)-dimcnsional submanifold of physical states consistent with the underlying internal symmetry. 

New constraints lead to decreasing number of available states for the system during its time evolution. The 
equilibrium distribution can be constructed by the Lagrange multipliers method. The multipliers related to the 
"central" constraints are well known chemical potentials. The multiplier related to the "nonabelian" constraint 
( [l^ ) is more complicated. Because the constraint is a nonlinear one, the corresponding multiplier cannot be 
treated as a standard additive thermodynamical potential. 
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